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Рассматривается алгоритмическая и программная реализация моделирования распределенных систем техно- 
логических объектов. Предлаг ается технология разработки иерархической структуры классов, сог‚ ласованной 
с множественным наследованием в языке Си++, при написании такого рода программных средств. Описы- 
вается исследование данного программного средства при моделировании типовых имитационных МОДЕЛЕЙ. 
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Введение. В оперативном управлении процессами транспортировки и распределения веществ 
в химико-технологических объектах и —производствах важное место отводится вопросу компью- 
терного моделирования распределенных в пространстве и времени технологических объектов. 
Необходимо использовать проверенные методы, обеспечивающие высокую точность моделирова- 
ния. Такие методы уже реализованы и развиты в пакетах МаЧаБ, Оскауе, 5сЙаЬ на основе библио- 
теки ОБЕРАСК. Вычисления в этих библиотеках поддержано программами матрично-векторных 
операций из библиотек (Ваз, [араск, Мпраск). 

В качестве примера имитационного моделирования сложных динамических систем с ие- 

рархической структурой средствами объектно-ориентирован-ного программирования взят участок 
магистрального газопровода. Имитационное моделирование подразумевает сохранение структуры, 
связи между компонентами, способа передачи информации [1], присущих моделируемой системе. 
Элементы системы заменяются их моделями, которые «обвязываются» информационными связя- 
ми так же, как и объекты в моделируемой системе. 
Область исследований. Компьютерное моделирование основано на математическом моделиро- 
вании и использует в качестве основы математические модели, тем не менее, математические 
предположения ослаблены для того, чтобы пользователь мог взаимодействовать с моделью 
интерактивно. Изучение модели пользователем в ходе интерактивного взаимодействия называет- 
ся имитацией. Имитация как метод решения нетривиальных задач получила начальное развитие в 
связи с созданием ЭВМ в 1950-1960-х гг. 

Выделяют две разновидности имитации: 

— метод Монте-Карло (метод статистических испытаний); 

— метод имитационного моделирования (статистическое моделирование). 

Цель имитационного моделирования — воспроизведение поведения исследуемой системы 
на основе анализа наиболее существенных взаимосвязей между ее элементами или, другими сло- 
вами, разработка симулятора (5/тшаНоп тоде/тд - англ. термин) исследуемой предметной об- 
ласти для проведения различных экспериментов. 

Имитационное моделирование - это имитация (симуляция) поведения системы во време- 
ни, которым можно управлять: замедлять в случае с быстропротекающими процессами и ускорять 
для моделирования систем с медленной изменчивостью. Можно имитировать поведение тех объ- 
ектов, реальные эксперименты с которыми дороги, невозможны или опасны. 

Математическая модель выражает существенные черты объекта или процесса языком 
уравнений и других математических средств. Модель учитывает возможности исполнительного 
устройства, получающего результаты из исходных данных, используя построенную модель и ал- 
горитм (рис.1). Математические формулы, позволяющие вычислить поведение модели, называют- 
ся вычислительной математической моделью. 


500 


Вестник ДГТУ, 2010. Т. 10. №4(47) 

















Огрубление объекта] 
(процесса) 


Определение целей 
моделирования 


Исюодный 
объект (процесс) 





Поиск математического 
описания 













Выбор метода 
исследования 













Анализ 
результатов 





Разработка алгоритма и 
программы для ЭВМ 





Отладка и тестирование 
програмы 





Рис.1. Компьютерное имитационное моделирование 


Компьютерное моделирование — реализация вычислительных процессов в ЭВМ на 
основе математических моделей для исследования реального процесса на основе вычисленной 
информации. 

Наиболее полный и однозначный переход от процесса к модели получается при воспроиз- 
ведении физических законов. Более сложные процессы, например, технологические, с опреде- 
ленной долей точности можно представить как суперпозицию этих законов. Для сложных процес- 
сов установить структуру физических законов, энергетических и информационных связей позво- 
ляет последовательный анализ взаимосвязей и принятие обоснованных допущений при абстраги- 
ровании. Поэтапное абстрагирование и рассмотрение «Переменных порядка» [2] системной ди- 
намики на разных уровнях абстракции приводит к построению иерархии математических моделей. 
В программировании такие переменные - ячейки памяти, через которые передается входная 
и выходная информация между имитационными моделями разных уровней абстракций, например, 
между реализацией метода решения дифференциальных уравнений и визуализатором состояния 
пространственно-распределенного объекта. 

При абстрагировании некоторые различия моделей отбрасываются, их объединяют в один 

класс, в рамках которого они не различимы, при этом некоторые важные свойства объектов те- 
ряются. Чтобы иметь возможность промоделировать процессы на объектно-ориентированном 
языке программирования во всей их сложности, надо, последовательно абстрагируясь от несуще- 
ственных свойств (на определенном этапе), строить иерархию классов, учитывающую отброшен- 
ные свойства объектов. Множественное наследование позволяет в результате получить объекты 
с нужными свойствами и необходимой степенью абстракции (синтез производных классов). Таким 
образом, задача создания среды объектно-ориентированного моделирования сложной системы 
включает этапы построения иерархии классов, порождения (придания объектам новых свойств) 
и их увязывания (формирования известных информационных связей). 
Результаты исследований. Система уровней абстракций помогает упорядочить процесс по- 
строения иерархии взаимодействующих объектов в парадигме объектно-ориентированного про- 
граммирования. В ходе анализа структуры магистрального газопровода было выделено шесть 
уровней, на которых особым образом выделяются классы объектов. Такое дробление на уровни 
подходит для всех распределенных технологических объектов и для их моделей, например, уча- 
стков магистрального газопровода, теплотрасс, химических реакторов, тепловых технологических 
печей, транспортеров сыпучих материалов. Для таких моделей технологических объектов можно 
описать все шесть уровней, а для компьютерной модели идеального вытеснения и идеального 
смешения - только первые четыре. Описание уровней приведено ниже. 
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Кибернетический уровень. Кибернетика - наука об управлении, связи и переработке ин- 
формации. Основным объектом исследования в кибернетике являются кибернетические системы 


(рис.2). 





Соп=Еаг 3 





Сопз ап 





Сопеаг 5 





Соп&ап 8 бибзузе тб 


Рис.2. Кибернетическая система на примере ячеечной модели 


Состояние элемента кибернетической системы может меняться как самопроизвольно, так 
и под воздействием тех или иных входных сигналов, получаемых им извне (из-за пределов рас- 
сматриваемой системы) либо от других элементов системы. В свою очередь, каждый элемент сис- 
темы может формировать выходные сигналы, зависящие в общем случае от состояния элемента 
и воспринимаемых им в рассматриваемый момент времени входных сигналов. Эти сигналы либо 
передаются на другие элементы системы, при этом служат для них входными сигналами, либо 
входят в качестве составной части в передаваемые за пределы системы выходные сигналы всей 
системы в целом. 

Ввод статических параметров осуществляется с помощью запроса у пользователя 
(интерактивно). В процессе работы происходит моделирование динамических процессов с одно- 
временным опросом реакции пользователя (см. рис.3). На выходе получаем визуализацию полу- 
ченных статических состояний. 
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Рис.3. Алгоритм интерактивной визуализации 


Механический уровень. В теоретической механике используется понятие фазового про- 
странства — переменные состояния (положения, скорости, ускорения материальных тел). В связи 
с этим в моделях нужно учитывать в несколько раз больше переменных, нежели количество на- 
копителей энергии в объекте моделирования. 

Абстрактные объекты, инкапсулирующие (объединяющие и скрывающие) все эти пере- 
менные и формирующие динамику их изменения, относятся к классу «Фазовое пространство». 
Причем размерность фазового пространства, объектов, определяемых в форме Коши, в несколько 
раз больше, чем самих переменных состояния, поэтому класс «Переменные состояния» является 
наследником класса «Фазовое пространство», учитывающего общие свойства динамических мо- 
делей. 

Физический уровень. Физическим этот уровень назван потому, что в его основу легла идея 
разделения классов, описывающих различные модели. В физике различают идеализации процес- 
сов — динамические, сосредоточенные в пространстве, распределенные в пространстве и стацио- 
нарные во времени. В связи с такой классификацией при моделировании сложных технологиче- 
ских объектов выделяют три вида абстрактных моделей. Объекты первого вида моделируют про- 
странственную сетку с постоянным или переменным шагом, необходимую для решения задач чис- 
ленных методов. Такие модели позволяют промоделировать пространственные размерности рас- 
пределенных объектов. Второй вид объектов моделирует суперпозицию нескольких воздействий, 
подаваемых в одну точку пространства, и распределение их по нескольким каналам. Такие без- 
размерные модели характеризуются параметрами — количеством входов и выходов. И третий вид 
— абстрактная модель, включающая алгоритмы учета хода виртуального модельного времени. 
Алгоритмы, реализуемые в этом классе, позволяют программам численного моделирования ре- 
шать задачи только на уровне численных методов и не заботиться о вычислении таких обобщаю- 
щих показателей, как процент выполнения вычислительной работы по моделированию и обнов- 
лению пользовательского интерфейса. Композиция этих видов порождает объекты полезные при 
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моделировании физических систем общей природы, без привязки к техническим особенностям 
или технологии протекания того или иного процесса. 

Математический уровень. Представляет алгоритмы аппроксимации непрерывного распре- 
деления параметров сеточными (табличными) функциями. При моделировании технологических 
объектов происходит замена моделей, заданных в непрерывном времени, дискретными моделями, 
реализуемыми компьютерными програм- 
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моделирования химико-технологического объекта 


называется  пространственнои —сеткои. идеального смешения 


Сетки бывают равномерные и неравно- 

мерные. Для построения неравномерных сеток, как правило, необходимо серьезное обоснование, 
так как на неравномерной сетке формулы аппроксимации усложняются. Ограничимся пока рас- 
смотрением равномерной по пространству сетки узлов. 

Пространственно-временная аппроксимация в узлах этой сетки, в частности, метод пря- 
мых или конечно-разностная аппроксимация пространственных градиентов, позволяет получить 
аппроксимированную математическую модель, описываемую системой дифференциальных урав- 
нений, которые можно решить с помощью различных численных методов. Наиболее частое при- 
менение находят методы Рунге-Кутта, Розенброка, Милна и т.д. К сожалению, не найдены науч- 
ные публикации, анализирующие совместимость этих методов с аппроксимированной на про- 
странственной сетке моделью распределенной среды. Соответствие исходной и аппроксимиро- 
ванной моделей достигается соблюдением дифференциальных законов протекания моделируемых 
процессов, а иногда интегральными зависимостями, описываемыми математической физикой. 
Между исходной и аппроксимированной моделью должна существовать тождественность, доказы- 
ваемая предельным переходом к бесконечно малым величинам по пространству и времени. По- 
строение компьютерной вычислительной модели происходит абстрагированием от непрерывной 
распределенности среды в конкретном аппарате. В общем случае принимается гипотеза о кусоч- 
ной непрерывности и-й степени. Также производится переход от непрерывного времени к дис- 
кретному, ход которого определяется конкретным численным алгоритмом и влияет на качество 
моделирования. Для качественной визуализации процессов в модели необходим обратный пере- 
ход от табличных функций, определенных на «дискретном времени», к непрерывному времени 
и пространству с точностью до компьютерной погрешности представления чисел с плавающей 
точкой. Такой переход реализуется межузловым восполнением — полиномиальной интерполяцией. 

Технический уровень. Описывает технические конструкции, являющиеся составными час- 
тями технологических объектов, т.е. виртуальные линейные участки (участки газопровода без 
поворотов), виртуальные сопротивления (сопротивления газопровода, связанные с изгибом тру- 
бы, заслонки, газоперекачивающие агрегаты). 

Технологический уровень. На технологическом уровне удобно объединить различные ви- 
ды технологических объектов в классы: участки газопровода, компрессорные станции. На данном 
уровне не существенны такие низкоуровневые аспекты реализации модели, как способ аппрокси- 
мации или шаг пространственной сетки. Разработаны различные программные модели (классы 
объектов), в некоторых случаях допускающие программную стыковку через интерфейсы 
приОшриё. Объекты классов моделируют: 

— СазР!рете - участки газопровода; 

— Сотргезвог${аНоп — компрессорные станции; 

— Неа{п9Мат — теплотрассы; 

— Сопуеуог — транспортеры сыпучих материалов. 

Моделирование при описанном подходе является результатом взаимодействия множества 
внутренних компонентов программы. На схеме (рис.5) можно увидеть четкую иерархию классов 
в соответствии с предложенными выше уровнями. 
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Рис.5. Иерархия классов: 
КС - компрессорная станция; РЗ -— запорно-регулирующая арматура; 
УМГП - участок магистрального газопровода; ДУЧП - система дифференци- 
альных уравнений в частных производных; ГПА - газоперекачивающий 
агрегат; МТМОбу$ет -— система со множеством входов и выходов 


Основной поток информации при моделировании циркулирует между компонентом Ргодесё, 
реализующим численный метод решения систем дифференциальных уравнений, относящимся 
к классу бо№ег, и набором объектов МодЕе|, отвечающих за вычисления правых частей дифферен- 
циальных уравнений (ДУ). Эти объекты взаимодействуют между собой при помощи интерфейсов, 
описанных в классе приОш ри. 

Такая схема взаимодействия мотивирована модульным принципом и позволяет подклю- 

чать внешние модули без перестройки программной и алгоритмической структуры разработанно- 
го программного модуля. При традиционном моделировании в среде пакетов инженерно-научных 
расчетов МаЧаб или СМУ Осбауе правые части системы ДУ вычисляются в пользовательской 
функции. Это обосновано, когда уравнения относительно простые, но когда расчет уравнений по 
вычислительной сложности сопоставим с работой численного метода, имеет смысл использовать 
специальные архитектурные решения. Предложенный подход значительно уменьшает вычисли- 
тельные затраты при моделировании объектов, поточно взаимодействующих непосредственно 
через переменные состояния. 
Применение результатов исследований в программировании имитационной модели. 
Разработанная система уровней позволяет стоить иерархию классов программных объектов, 
а затем связывать объекты в ходе выполнения программы во время моделирования. Структура 
«моделирующего дерева» в простом случае вырождается в список взаимодействующих моделей 
разных уровней абстракции, например, структура имитационной модели участка магистрального 
газопровода (рис.6). Предложенный подход к декомпозиции программ и алгоритмов на модули 
позволяет упорядочить разработку программ, моделирующих сложные технологические процес- 
сы. 
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Рис.6. Взаимодействие программных объектов при имитационном моделировании 


обобщенный алгоритм компьютерного моделирования распределенных технологических 
объектов набором взаимодействующих имитационных моделей и вспомогательных модулей со- 
стоит из следующих шагов: 

1. Определяются параметры, необходимые для работы каждой из имитационных моделей 
и вспомогательных модулей. Данные вводятся пользователем-исследователем или извлекаются 
из базы данных. В частности определяются параметры уравнений, которые явно не зависят ни от 
времени, ни от пространственной координаты. 

2. Вспомогательные алгоритмы в модуле Арргохитаюог определяют размещение простран- 
ственных узлов для каждого из моделируемых аппаратов и их внутреннего пространства. Выде- 
ляется память под переменные состояния и фазовые координаты. 

3. Инициализируются параметры численного метода, выбранного пользователем. Для па- 
кетного (неинтерактивного) моделирования задается наибольший Л„„„, наименьший й„„ шаг ме- 
тода И [нь (нах —- начало и конец временного отрезка, на котором интегрируется уравнение, вели- 
чина локальной точности метода =. 

4. Для непрерывного интерактивного моделирования задается допустимая точность вы- 
полнения (невязка) интегральных балансных уравнений. 

5. По уравнениям статики рассчитывается комплексный режим статики всех имитацион- 
ных моделей, входящих в систему. 

6. Дальнейшие вычисления проводятся в цикле. Итерации выполняются для каждого дис- 
кретного момента времени, определяемого численным методом. Для вычисления скоростей изме- 
нения признака (концентрация, плотность, масса, скорость) во времени применяется объект Ар- 
ргохитаюг, вычисляющий приближенно пространственные градиенты признаков и использующий 
модель конкретной сплошной среды, заданной дифференциальным уравнением в частных произ- 
водных (ДУЧП). Пространственные градиенты используются в последней модели для вычисления 
скоростей изменения признака во времени. В свою очередь, эта информация возвращается мето- 
ду интегрирования по времени. При расчете параметров в граничных ячейках задействуется ин- 
терфейс, определенный в классе приОшфри{, для согласования граничных условий с соседними 
технологическими агрегатами. 

7. Найденные скорости изменения признака в численном методе обобщаются в «Вектор 
наклона интегральных кривых», /. Программная реализация численного метода (Милна, Розен- 
брока, Эйлера, Рунге-Кутта) вычисляет новую точку состояния системы, координаты которой со- 
храняются в реализации класса Зайебрасе. 

8. Цикл повторяется до тех пор, пока не достигнет конца временного отрезка ({„„„. Выход 
из цикла производится в случае выхода модели в критический режим, в пределе называемый син- 
гулярностью. Когда производится непрерывное интерактивное моделирование (симуляция), поль- 
зователь-исследователь имеет возможность приостановить процесс моделирования или завер- 
шить его, сохранив хронологию изменения переменных состояния в базе данных (см. рис.3). 
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